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ABSTRACT 


This  paper  is  concerned  with  (i)  the  maximization 


of  net  return  £  n.(y  ,t  )  and,  (ii)  the  maximi- 
j-1  J  J  J  n 

zation  of  average  net  return  £  tt  Cy,  ,t)  .  of  a 

j-1  j  J 

deterministic  multicommodity  system  subject  to 


linear  restrictions  on  the  inventory  levels  y^  and 
the  review  periods  t^  (J*’ ,  . . . ,  n)  . 

The  first  problem  is  a  symmetric  quadratic  program. 
For  the  second  problem,  an  algorithm  is  given  for 
finding  the  optimal  solution. 
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I.  INTRODUCTION 


1.  In  this  paper  Is  considered:  the  maximization  of  the  net  return 

J  »j(yj,tj)  and  the  average  net  return  ][  °*  *  deterministic 

multicommodity  system  subject  to  linear  restrictions  on  the  inventory  levels 
y^  and  the  review  periods  t^  . 

2.  It  will  be  useful  to  give  first,  a  short  description  of  the  deter¬ 
ministic  multiconroodity  system  whose  the  net  return  is  n(y,t)  .  The 
characteristics  of  such  a  system  are  the  following  [7]. 

(1)  Demand  is  deterministic  at  a  constant  rate  of  R  quantity  units 
per  time  unit,  R  >  0  . 

The  replenishment  rate  is  infinite. 

$  . 

h  -  inventory  holding  cost,  h  >  0,  dimension:  unit‘  Unit- time 
c  ■  cost  of  production,  dimension:  '^j^r 
r  •  revenue,  dimension:  $ 

K  *  ordering  cost,  dimension:  $ 


(ii) 

(iii) 

(iv) 

(v) 
(Vi) 
(vii) 


P  -  penalty  (shortage  cost),  dimension:  Unit  .Unit_Sho'r't 


(viii)  P  -  penalty  (good  will  lost),  dimension:  -Short 

(ix)  y  ■  inventory  level,  dimension:  Unit. 

(x)  t  ■  review  period  time,  dimension:  Time 

(xi)  ir(y,t)  -  net  revenue,  dimension:  $ 

(xil)  h,  c,  r,  k,  P  and  P  are  nonnegative  constans,  whereas  y  and  t 

are  variables. 

The  value  of  the  net  return  w(y,t)  is  given  by 


(1) 


*(y.t) 


F(y,t) 

C(y,t) 

H(y,t) 


if  y  <  0 

if  0  <  y  <  Rt 

if  y  >  Rt 
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where : 

(2)  F(y,t)  -  (r-c)Rt 

(3)  G(y ,t)  -  (r-c)Rt 

(4)  H(y,t)  -  (r-c)Rt 


[K  +  p(fe  ~  y)  +  pRt] 

[K  + fey2  +  fe  (Rt_y)2  +  p(Rt-y) 

[k  +  h(y  -  f  t)tj 


From  (2),  (3)  and  (4)  it  can  be  observed  that  ir(y,t)  has  the  following 
properties : 

(a)  :< (y , t )  is  a  continuous  function  of  y  (unrestricted)  and  t  >  0 

(b)  exists  and  is  f.  continuous  function  except  at  the  points 
ay 

(0,t)  and  (Rt , t)  . 

(c)  '^3^ ^  exists  and  is  a  continuous  function  except  at  the  points 

(Rt,t)  . 

(d)  n (y , t)  is  a  concave  function  for  all  y  and  t  such  that 

0  <  y  <  Rt  . 

Define  the  negative  net  return  N(y,t) 


(5) 

N(y,t)  -  -ir(y,t) 

(6) 

N(y|t)  -  N(y,t  - 

From  (2),  (3)  and  (4) 


pt  i  y  <  0 

lN-(y.l.O.  .  i  y  -  I  (Rt-y)  -  p  if  0  <  y  <  Rt 
7 


3N(y|t) 

3y 


y-Rt- 


if  y  >  Rt 
ht-p 
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The  graphical  representation  of  N(y|t)  is  given  in  Figure  1. 


It  is  clear  that  the  function  N (y | t )  attains  its  minimum  at  a  point  y* 

*  p  *  p 

such  that  y  «  Rt  if  t  <  —  or  0  <  y  <  Rt  if  t  >  —  • 

h  -  J  -  h 

Lost  Sale  Case 

In  the  model  discussed  before  it  was  assumed  that  all  demands  incurred 
when  the  system  was  out  of  stock  were  backordered.  If  demands  occurring  when 
the  system  is  out  of  stock  are  lost  forever,  the  value  of  the  net  return 
"(y,t)  will  depend  on  the  length  of  time  for  which  the  system  is  out  of  stock. 
In  this  case,  since  backorders  are  not  allowed,  the  penalty  due  to  shortage 
is  null,  whereas  the  penalty  due  to  the  good  will  lost  is  equal  to  r-c,  i.e., 

P  *  0  and  p  -  r-c  [9];  with  these  changes  in  the  equation  (1),  the  net 
revenue  for  the  Lost  Sale  Case  becomes 

0  if  y  <  0 

(r-c)y  -  [k  +  y2J  if  0  <  y  <  Rt 

(r-c)Rt  - [ K  +  h(y  -  ~  tjt]  if  y  >  Rt 


(8) 


<r(y,t)  - 
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II.  RESTRICTIONS 

1.  Restrictions  on  the  Inventory  Levels 

There  can  be  many  sorts  of  interactions  between  items  when  an  inventory 
system  stocks  many  items.  For  example,  warehouse  capacity  may  be  limited,  in 
which  case  the  items  compete  for  the  floorspace;  or  there  maybe  an  upper  limit 
to  the  maximum  investment  in  inventory  whence  the  items  compete  for  investing 
dollars;  or  the  total  weight  of  items  that  the  warehouse  floor  can  tolerate 
may  be  limited;  etc. 

Consider  the  case  where  there  is  an  upper  limit  b^  to  the  square  feet 
of  warehouse  floor  space.  Suppose  that  n  items  are  being  stocked  and  one 
unit  of  item  j  takes  up  blj  square  feet  space,  then  if  the  floor  space 
constraint  is  not  to  be  violated  at  any  time,  it  must  be  true  that 
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In  some  models  all  demands  must  be  met  from  Inventory  so  that  no  back¬ 
orders  or  lost  sales  are  allowed,  in  this  case  the  restrictions  are 

yi  m  Yj  J-1 . .  l7l# 

If  the  maximization  of  £  *J (Yj ’tj>  i®  on^y  f°r  one  period,  then  all 


J-l 


shortages  are  lost  sales  and  the  optimal  solution  must  be  such  that 

*  *  *  * 

y^  <  Rjtj  j-1,  •••»  n;  that  Is,  y^  >  R^t^  must  not  permitted  since  there 

*  * 

Is  not  a  second  period  for  selling  the  difference  y^  -  R^t^  .  Therefore, 
y,  <  R.t,  must  be  s  restriction  for  all  Index  1-1,  ...,  n  . 

J  -  J  J  J 


6 


Program  A' 


7 


(11) 


n 

Minimize  N(y,t)  «  [  N  (y  ,t  ) 

1=1  J  J  J 

Subject  to  b  y  +  . .  .  +  b,  y  <  fa¬ 
ll  1  In  n  ■  1 


b  ,y,  +  . .  .  +  b  y  <b 
rl7l  m7n  ■  r 


d-~t.  •••  d-  t  <  d. 

11  1  In  n  ■  1 


0 :  yj  ;  Vr  J‘1 . n  • 

where 

y  “  y(yx.  •••.  yn)',  t  -  (t1,  tn)',  (y,t)  =  (y1 . y^;  tj,  ....  t  ) 

all  coefficients  of  y^  and  t  ,  and  b^  and  are  constants  for  all 

J-l,  n;  i-1,  ....  r;  k*l,  s  .  6(y^,0)  is  defined  by  (9)  and 

bij  -  °»  bi  >  0  for  a11  i*’1»  •••»  r;  J"1 . n  • 

Theorem  1 

*  * 

The  solution  (y  ,t  )  is  optimal  for  the  program  A  if  and  only  if  it 
is  optimal  for  the  program  A'. 

Proof : 

Define 

N  ®  ( 1 ,  •  •  *  i  n}  • 

S  -  {(y,t)|(y,t)  satisfies  (10)}  . 
s'  ■  {(y»t) ! (y,t)  satisfies  (11)}  . 

Let  (y°,t°)  be  a  solution  of  (10),  i.e.,  (y°,t°)c  S.  Associated  with  (y°,t°) 


dcf ined  the  sets  : 
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Nj  -  (j I jcN,  y° 

<  0) 

• 

N2  -  (j|jcN,  0 

0 

-  yJ 

-V 

n3  "  (J IjcN,  y° 

o  •  ’ 

(y  *  t  )  “  (y^  »  • • • »  yn; 

tl»  " 

'  '  »  ' 

c°)’ 

n 

6uch  that 

0 

if 

jeN1 

»i-j 

yj 

if 

jeN2 

| 

i  Yj 

if 

jeN3 

The  vector  (y',t°)  satisfies  the  restriction  (11),  because  for  all 

\ 

i“l . .  it  is  true  that 


biJyj  ’  biJ 

6(y°,0)  - 

0  if 

jeNj 

biJyj  "  bU 

o 

yj 

if 

JeN2 

Vj  -  bij 

0 

yJ 

if 

JcN3 

Furthermore,  since  yj  >  y° 

if  jcN^ 

’i-'- 

lf 

jeN2  and 

yi  "  yj  if 

JeN^,  a  direct  application  of  (7)  gives 

• 

V’i-1? 

i N:  <yrt°) 

if 

lzU1 

Y'r'? 

■ 

if 

jcN2 

VY'? 

if 

jeN3 

That  is,  for  each  solution 

(y0,t°)c  S 

there  is  a 

solution 

(y ' ,t°)e  S' 

such  that  N(y',t°)  <  N(y°,t°)  .  This  implies: 


minimum  M(y,t)  <  minimum  N(y,t) 
(y,t)e  S'  (y,t)e  s 


(12) 
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On  the  other  hand 

S'  C  S,  implies: 


(13) 


minimum  N(y,t)  >  minimum  N(y  t) 
(y.t)e  s'  “ 


Fr°m  (12)  and  (13),  the  Theorem  follows. 


Q.E.D. 


In  view  of  Theorem  1.,  it  will  be  sufficient  to  consider  the  function 

N(y,t)  only  on  the  region  0  <  y  <  r  t  4.1 

-  yi  c  J  j*  J  1 . .  • 


Then,  from  (3)  when  0<y  <Rt  K  - 

-  i  -  J  J  ’  j  0 

n 


N(y,t)-  l  N  (y  t) 
J-l  J  J  J 


(14) 


n 

’.I 

J-l 


h-  2  P 
+ 


'j  T  2R~  (Vi~Vir 


^iVyVVrYj’ 


using  matrix  notation,  the  first  te™  on  the  right  hand  side  can  be  written 
in  the  form 


(15) 


where 


n 

I 

J-l 


H-  2  ^ 
y 4  + 


VJ  ^  *r  (V)'yir 


-  (y.O'Q(y,t) 


(y,t)  -  (yx,  •••.  yn;  tx>  ...,  tQ) '  and  Q  is  the  2n  x  2n  matrix 
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Vpi 

2R, 


Q  - 


h  +P 
n  n 

2R 


R  P 
n  n 


Observe  that  the  matrix  Q  is  positive  (semi-)  definite  since  all  coefficients 
2  2 

of  yj  and  (R^t^-y^)  on  the  han<*  side  of  (15)  are  nonnegative.  That 

Is,  (y,t)'Q(y,t)  >  0  for  each  vector  (y,t).  In  particular,  the  left  hand 


side 

of  (15) 

is  always  positive  for 

each  vector  (y,t)  i  0,  if 

0 

for 

•11 

J-l,  .. 

. ,  n;  i.e.,  the  matrix 

Q 

is  positive  definite  if 

* 

0 

for 

•11 

J-l,  .. 

•  i  n  • 

(Note 

:  at  the  beginning  it  was 

assumed  h^  >  0,  R^  >0 

and 

>  0 

m 

for  all  J-l,  ....  n) . 

Then,  an  equivalent  formulation  of  the  program  A'  is 


(16) 


where 


Minimize  N(x)  -  x’Ax  +  g'x 
Subject  to  Ax  +  b  >  0 

x  >  0 

* 


x  -  (y,t)  *  (y^ ,  ....  yn;  ^ . tfl) 


g  «  (-Pi »  •••.  ”*_»  (Pi+c.-r. )R.  ,  ....  (P  +c  -r  )R  )' 
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and  b  is  the  (r+S+n)  column  vector 

b  •  (bj,  ••••  b^j  d^,  ...,  d^;  0,  ...,  0) 

2.  Existence  of  the  Solutions. 

The  set  S'  -  {(y,t)|(y,t)  satisfies  (11)}  is  a  closed  convex  set. 

If  in  addition  it  is  bounded,  then  the  program  A'  always  has  an  optimal 
solution,  because  N(y,t)  is  a  continuous  function  and  a  continuous  function 
defined  on  a  compact  (closed  and  bounded)  set  has  an  absolute  minimum  on  It. 
Let  and  be  such  that 

S2  *  {y| (y,t)e  S')  and  S2  •  (t|(y,t)c  S'} 

The  set  is  clearly  bounded  because  in  the  set  of  restrictions  (11)  it 

was  assumed  tnat  >  0  and  b^  >  0  for  all  1*1,  . ..,  r;  J»l,  . ..,  n  . 

The  set  may  be  unbounded. 

However,  if  >  0  for  all  J»l,  ...,  n;  then  the  program  A'  has  always 
an  optimal  solution.  This  affirmation  follows  immediately  from  the  following: 
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(i)  S  is  a  closed  convex  set,  and  (ii)  when  >  0  for  all  j«l,  . ..,  n 

the  objective  function  N(y,t)  cannot  be  unbounded  in  the  direction  of  ex- 
tremization,  since  it  is  the  sum  of  a  positive  definite  quadratic  form  and  a 
linear  expression 

3.  There  are  many  ways  t'  obtain  the  optimal  solution  of  program  A',  a 
quadratic  symmetric  program.  One  method  is  the  following:  See  [1], 

Find  Z  >  0,  such  that 

m 

MZ  +  q  >0  and 
Z'(Mz  +  q)  -  0 

where 

M  «  ^2Q  -  A 'y,  q  -  (g,b)  and 

Z  •  (x,X)  .  The  components  of  the  vector  X  are  the  nonnegative  multipliers 
associated  with  the  inequalities  of  (16). 


Example  1 


A  retail  merchant  in  city  A  makes  a  weekly  trip  to  city  B  in  order  to 
refill  his  supplies.  The  truck  he  uses,  restricts  the  volume  of  the  goods. 
Therefore  the  merchant  must  decide  how  much  of  each  commodity  he  should  take 
such  that,  (i)  the  truck  restrictions  are  not  violated,  and,  (11)  his 
total  profit  is  maximized. 

All  shortages  are  lost  sales  and  the  set  up  cost  is  the  cost  of  the  trip 


Cj  is  the  buying  price  at  city  B;  r ^ ,  h^  and  are  respectively  the 

selling  price,  inventory  holding  cost  and  demand  of  iters  j»l . .  v 


and  w^  are  the  volume  and  weight  of  item  j  .  Finally  V  and  W  are  the 
volume  and  weight  capacities  of  the  truck. 


For  instance,  with  n=2  and  the  following  values: 


Item 

rJ 

CJ 

hj 

PJ 

Vj 

WJ 

K 

V 

U 

1 

5 

3 

0.4 

20 

0 

2 

6 

3 

100 

2400 

1500 

2 

3 

2 

0.5 

50 

0 

1 

4 

5 

For  reason  of  space  all  dimensions  were  dropped. 

The  unit  time  is  1  day, 

"  "  volume  is  1  cubic  foot, 

"  "  weight  is  1  pound. 

Then,  the  problem  is 

Minimize  [ife  +  200  '  2yl  “  y2  *  100 

Subject  to  3y^  +  5y2  <  1500 

6y^  +  4y2  <  2400 
y1  <  20x7 
y2  <  50x7 


The  optimal  solution  is 


100,  y2  -  100 
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n 


IV.  MAXIMIZATION  OF  £  "4(y.,t)/f. 

j-1  J  J 

The  maximization  of  the  average  total  net  revenue  per  unit  time  when  pll 

replenishments  are  made  up  simultaneously  and  the  inventory  levels  are  subject 

to  linear  restrictions  is  considered  in  this  section. 

n 


For  convenience,  one  will  minimize  £  N.(y.,t). 

j-1  3  3 


J-1 


jvyJ’w/t  ‘ 


Instead  of  maximizing 


1.  The  Program  Under  Consideration  Is  Program  B: 

n 

Minimize  N(y,t)^-  [  ^^j’^/t 

Subject  to  b^  6(y^,0)  +  ...  +  bj  6(yn>0)  <  b ^ 


brl6(y1,0)  +  ...+brn«(yn.0)  <br 

y^  is  unrestricted  J-1,  ...,  n 

t  >  0  . 

m 

vhe  re 

(y,t)  -  (y^,  ...,  y^.t)'  Is  an  (N+l)  column  vector  and  the  restriction  on 
the  inventory  levels  y^  are  the  same  as  in  program  A*. 


Consider  also 


15 


Program  B' : 


Minimize  N(y,t)^fc-  \  Nj  »t) /t 
Subject  to  +  ...  +  blnyn  <  b][ 


b  ,y.  +  . . .  +  b  y  <  b 
rl7l  nrn  ■  r 


0  £  ^  J-1 . n 

t  >  0  . 


Theorem  2 


*  * 


The  solution  (y  ,t  )  is  optimal  for  the  program  B  if  and  only  if  it  is 
optimal  of  the  program  B'. 

With  tj«=t  (j«lf  ....  n)  the  proof  of  this  theorem  follows  the  same 
line  as  the  one  for  Theorem  1. 

Define 


(19) 


T  *  ((y,t)|(y,t)  satisfies  (18)}  . 


Theorem  3 


*  * 


If  )  <  o  and  N(y+,t+)  <  0  then  t  <  t+  . 

t 

Where  (y  ,t  )  is  the  optimal  solution  of  program  Bf  and  N(y  ft  ) 

minimum  N(y,t)  . 

(y.t)e  T 


Proof : 


*  * 


If  (y  ,t  )  is  an  optimal  solution  of  program  B',  then 

»<y*.t*)  .  »(Y.t) 

t*  -  ‘ 


for  all  point  (y,t)c  T, 
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in  particular: 

do)  ,lnc.  (y+>t+)t  T  . 

t+ 


*  * 

r_ 

t 

On  the  other  hand, 

+  .+ 


N(y  ,t  )  <  N(y,t)  for  all  point  (y,t)c  T, 


in  particular: 

(21)  N(y+,t+)  <  N(y  ,t  ),  since  (y*,t*)e  T  . 


Using  the  hypothesis  of  the  theorem  in  (20)  and  (21)  one  can  get,  respectively 

(22)  ^  ^  >  —  and  (23)  V  <  1 

N(y  ,t  )  "  t  N(y  ,t+)  " 

From  (22)  and  (23) 

&T  <  *  1  •  “'■>'«  t‘  <  t+  . 

t+  -  N(y\t+)  ' 

Q.E.D. 

The  economical  meaning  of  this  result  is:  if  the  optimal  total  net 
revenue  for  the  one  period  analysis  is  positive  and  its  optimal  review  period 
_is  t+,  when  the  optimal  average  total  net  revenue  per  unit  time  is  nonnegative, 
then  t+  is  an  upper  bound  for  the  optimal  review  period  of  the  Steady  State 
Case. 


In  view  of  Theorem  2.,  it  will  be  sufficient  to  consider  the  function 
N(y,t)/t  only  on  the  region  0  <  y^  <  R^t  (J»l,  . . . ,  n)  . 

From  (3). 


N(y-,)/t  ■  Jj[Kl  +  2R^  »)2  +  2^  W^Wrj’VVl 


/ 1 


N* (y  ,t)  +  N"(t) 
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where 

n 

(24) 

N'(y.t)  =  l 
j-1 

-(pj  +  TL) yjj 

n 

K .  R.P . 

(25) 

N"(t)  -  l 

J-1 

t  +  (P^-c.^ 

*  • 

Define 


(26) 


T(t)  -  {y|(y,t)e  T}  . 


Note  that  T(t)  Is  the  set  of  solutions  y  of  (18)  for  a  fixed  value  of  t, 
and  T(t)  Is  a  compact  set. 

Then,  an  equivalent  formulation  of  the  program  B'  is 

Minimize  ffXyi*? 

Subject  to  (y,t)c  T  . 


Observe  that 


*  (27) 


Minimum  ^ 

(y,t)e  T 


minimum  / minimum  ^  j 

t  >  0  \y  e  T(t)  C  / 

minimum  /minimum  (N'(y,t)  +  N"(t)) 
t  >  0  lye  T(t) 

minimum  jN"(t)  +  minimum  N'(y,t)j 
t  >  0  \  ye  T(t)  / 


) 


Define 


(28) 

H(t)  *  minimum 

N  (y.i  0- 

y  €  T(t) 

u 

(29) 

h(t)  *  minimum 

N'(y.t) 

y  c  T(t) 

h(t)  +  N"(t) 
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Theorem  4 

The  program  B'  always  has  an  optimal  solution. 

In  order  to  prove  Theorem  4.,  the  following  definitions  and  theorems  are 
necessary. 

Definition  1.  Scalar  function 

A  correspondence  which  assigns  a  single  scalar  (a  real  number,  a  point 
in  R')  to  each  point  x  of  a  set  T  is  called  a  scalar  function  of  x  and 
is  written  as  0(x)  .  The  scalar  function  0(x)  is  said  to  be  defined  on  T  . 

Definition  2.  Ouasi-convex  function 

A  scalar  function  G(x)  defined  on  a  convex  set  TCRn  is  said  to  be 
quasi-convex  on  T  if  0(ax  +  (l-a)y)  <  max  (0(x),  0(y))  for  each  pair  of 
points  x,  y  e  T  and  all  real  a  >  0  . 

Definition  3.  Strictly  guasl-convex  function 

A  scalar  function  0(x)  defined  on  a  convex  set  FCRn  is  said  to  be 
strictly  quasi-convex  on  T  if  for  each  pai"  of  points  x,  y  e  T  such  that 
0(x)  t  0(y),  it  is  true  that  Q(ax  +  (l-a))<  max  (0(x) ,  0(y))  for  all  real 
a  >  0  . 

Theorem  5 

The  program:  minimize  N’(y,t)  such  "hat  y  e  T(t),  always  has  an  opti¬ 
mal  solution. 

Proof : 

Since  N'(y,t)  is  a  continuous  function  and  T(t)  is  a  compact  set, 
N'(y,t)  has  an  absolute  minimum  on  T(t)  . 
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Theorem  6. _ [3j 

The  function 


M(t)  «  minimum  N(y,t),  Is  convex 
y  c  T(t) 


Proof : 

The  function  N(y,t)  is  continuous  and  the  set  T(t)  is  compact. 
Therefore  for  each  number  t°  >  0  there  is  a  point  y°e  T(t°)  such  that 
N(y°,t°)  is  an  absolute  minimum. 

Let  t^  >  0  and  >  0,  with  t^  i  t2  •  Then  there  exists  y^e  T(t^) 


2 

and  y  e  T(t2) 

such 

that 

(30) 

M(t  ) 

-  minimum  N(y,t.)  -  N(y1,t.) 

X 

y  e  T(tx) 

(31) 

M(t,) 

2 

»  minimum  N(y,t-)  -  N(y  ,t.) 

y  c  t(t2)  1  1 

On  the  other  hand  the  set  of  restrictions  (19)  are  linear.  Therefore,  if 
y^  T(t1)  and  y2e  T(t2)  then  (ay1  +  (l-a)y?)e  T(at1  +  (l-a)t2)  for  all 
real  a  >  0  . 

Hence 

M(at1  +  (l-a)t2)  -  N(y ,at^  +  (l-a)t2> 

yeT(at j  +  (l-a)t2> 

<  N(ay1  l-a)y2,  a^  +  (l-a)t2> 

But  the  function  N(y,t)  is  convex,  then 

M(at^  +  (l-a)t2)  <  aN(y1,t1)  +  (l-a)N(y2 , t2>  . 


Using  (30)  and  (31) 
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Q.E.D. 


21 


Remark 

Since  the  function  M(t)  is  defined  and  convex  on  [0,  +*»)  it  is 
continuous  for  all  t  c  (0,  +“)  .  The  function  H(t)  *  is  defined  for 

all  t  >  0  and  it  is  the  ratio  of  two  continuous  functions,  hence,  H(t)  is 
continuous  function  for  all  t  >  0  . 


Proof  of  Theorem  h 


For  definition  (28) 


H(t)  *=  minimum 


N(y,t) 


y  c  T(t)  c 
1  n 

minimum  —  £ 

t  j-1 


K. 


* 


(R  t-y)+(P^+Cj-r  )R  t-P 


since 


Yj  “  RjC  *  *  *  ’ 


lim  H(t) 

t-»0+ 


n 

bounded  if  £  K  ■  0 

J-l  J 

n 

if  l  K  >  0  . 

j-l  5 


In  the  first  case  define  H(0)  •  lim  H(t)  . 

t-»0+ 

Then  the  function  H(t)  is  continuous  for  all  t  >  0.  In  Theorem  3, 

* 

it  was  showr  that  the  optimal  review  period  t  of  program  B  is  bounded 
by  t+  .  Hence,  since  H(t)  is  continuous  on  the  compact  set  [0,  t*] ,  the 
theorem  follows  immediately. 

If  H(t)-H-®,  one  can  find  a  number  {  >  0  such  that  H(t)  >  H(t+)  for 
t**0+  + 

t  <  £  .  But  H(t)  is  continuous  on  the  compact  set  [£,  t  )  .  Hence  H(t) 

has  an  absolute  minimum  on  t+]  . 


Q.E.D. 
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In  order  to  find  the  optimal  solution  of  program  B'  an  algorithm  is  given. 
The  algorithm  is  based  mainly  on  the  following: 

(i)  For  each  t  >  0,  the  program 

Minimize  N'(y,t) 

Subject  to  y  e  T(t) 


always  has  an  optimal  solution. 

(ii)  H(t)  is  a  continuous  (strictly)  quasi-convex  function  of  t  . 

3.  Algorithm 
Step  0 

(a)  Find  the  upper  bound  t+ 

(b)  Let  So  -  (0,  t+) 


Step  1 

Find  t.  such  that  H(t.)  ■  min  H(t) 

*  ^  ^  -  e 


Step  2 

(a)  Find  t'  and  t"  such  that 


t'  ■  min  [(t.-t)  > 
t£Si 


0] 


t" 


min  [(t-t.)  > 

tc^ 


0] 


if  tt  <  t+ 
if  t±  -  t+ 


(b)  Let  t't  -  j  (t,+t1)  and  t"1  -  \  (t'V^) 


Step  3 


Calculate  H(t'^)  and  H(t"^) 
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Setp  A 

Let  Si+1  "  SiUtt'i  ’  1  i* 

and  return  to  Step  1,  with  i+1  replacing  i  . 

In  Step  0,  the  upper  bound  t+  is  found  by  solving  the  program  (minimize 
N(y,t)  subject  to  (y,t)c  T) ,  see  Theorem  3. 

In  Step  3,  the  values  of  H(t'^)  and  H(t"^)  are  found  in  accordance 
with  the  formulas  (28)  and  (29). 

It  must  be  noted  that  at  each  iteration  it  is  necessary  to  solve  the 
following  programs 


Minimize  N'(y,t'^)  and  Minimize  N'ty.t"^) 

Subject  to  y  c  T^'^)  Subject  to  y  c  T(t"^)  . 

Define 


H(0)  -  lim  H(t) ,  and 

t-»0+ 

H(t  )  •  min  H(t) 
te[0,t  ] 


Theorem  7 


At  the  ^^(K-l,  2,  ...)  iteration  it  is  true  that 


(a)  t  -t'  -  t  and  t"-t  ■ 

IT  —  IT 


—if  t  <  t+ 
2*  K 


0  if  t  -  t 

K 


(W  . 


Proof : 

By  induction  over  K, 

(i)  the  theorem  is  clearly  true  for  K»1 
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(ii)  assume  that  it  is  true  for  K=r  and  consider  first,  the  case 
<  ‘+  • 

Then,  at  the  r  iteration 

+  .+ 


(32) 

where 


t  -t'  -  ,  t"-t  -  —  and  |t°-t_|  <  —  , 


r  2r 


r'  ■  2r 


(33) 


t'r  “  |  (t'+tr)  and  t"r  -  |  (t"+tr) 


Sr+1  "  SrU{t’r’  *V  * 


At  the  (r+l)**1  iteration 


Step  1 


H(t  ..)  -  min  H(t)  -  min(H(t  ),  H(t’  ),  H(t"  )) 
m  r  r  r 


(34)  Suppose  that  t^^  -  tf(i.e.  H(tr)  <  H(t’r),  h(t"r)) 


Step  2 
(35) 


t'  -  t'r,  t"  -  t"r 


Then,  from  (32),  (33),  (34),  and  (35), 

W  -  Vl't  ■  V  7  ■  7  <«*-*'>  *  7 ' 7  ■  • 

t"~tr+l  *  t"r”tr  '  7  '  7  *  7  '  £  ‘  • 


Hence,  the  part  (a)  of  the  theorem  is  true  for  all  K»l,  2,  ... 

Also,  t  e[t'r,  t"  ] ,  because  if  t  <  t'r,  then  there  exist  X  >  0, 
that  (Xt*  +  (1-X)tr)  ■=  t’r  and  H(Xt*  +  (1-X)tr)  -  H(t'r)  . 


such 
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But  the  quasi-convexity  of  H(t)  implies  that  HCt'^)  «■ 

H(Xt  +  (l-Xjt^.)  <  max(H(t  ),  H(tr>)  *■  H(tr>,  which  contradicts  the  assump- 

„  * 

tion  (34) .  Hence,  t  r  <  t 

In  a  similar  way  the  assumtion  tM  <  t  also  leads  to  a  contradiction. 
Hence  t  e[t'r,  t"^]  .  This  result  and  the  formulas  (32)  and  (33)  imply  that 


Hence,  | t  — t  <  1  .  Then  part  (b)  of  the  theorem  is  true  for  all 

’  ‘  r+1'  «=  2r+^ 

positive  integer  K-l,  2,  ... 

A  similar  analysis  for  the  cases  t  ^  -  t'  or  t"  in  (34),  completes 

the  proof  with  the  assumption  that  t  <  t+  . 

Following  the  same  steps,  the  analysis  of  the  case  t  ■  t+,  completes 

the  proof. 

Corollary  1 

lim  H(t  )  -  H(t*) 

The  corollary  follows  immediately  from  the  facts:  H(t)  is  a  continuous 

* 

function  and  11m  t  ■  t 
K-**  * 


Graphical  example 
,  H(t) 


Iteration  0 

(a)  Find  the  upper  bound  t+ 

(b)  S  -  {0,  t+} 

o 

Step  I 


H(t+)  ■  min  H(t) 
teS 

o 


Step  2 


(«)  t'  -  0,  t"  -  t+ 

(b)  ■  r  •  ■  t+ 

Step  3 


Calculate 


H(t+) 
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S  t  op  4 


si  ■  *.4r  •  t+ 


0  t+  +( 

u»  2  >  t  ( »  return  to  Step  1, 


Iteration  ? 


Step  1 


H(t  )  «  min  H(t)  «  H (—),  t  => 

frS  V  /  1  2 


(a)  t'  -  0,  t"  -  t+ 


Step  3 


Calculate  H ( L  ’  x )  and  HCt'^) 


Step  4 


s2  -  SjUft'j,  fp 


„  1  L  3t+  +1 

’4  ’2  ’  *  c  (•  return  to  Step  1. 


Iteration  112 


H(t2)  -  min  H(t)  -  H^-^,  c2  “  4“ 


Step  2 


(a)  t'  -  0,  t"  -  t  - 


a,  f2  -i  «+t2)  f2  ~  +  f).i£ 


Step  3 


Calculate  H(t'2)  and  H(t"2) 
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Step  4 


c 

w  3 


s2U{t'2,  t"2} 


return  to  Step  1. 


etc. 


* 
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